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I.  Introduction 

During  the  past  several  years,  the  area  of  laser* induced  surface 

chemistry  (LISC)  has  become  an  established  interdisciplinary  pursuit  among 

researchers  in  the  fields  of  physical  chemistry,  material  science  and 

optics.  LISC  has  been  reviewed  by  several  articles  for  various  applications 

1-3  3 

including  microelectronics  fabrication,  heterogeneous  catalysis  and 

4 

spectroscopy.  The  fundamentals  of  LISC  for  conventional  surfaces  were 

explored  in  our  earlier  review  papers, and  more  recently  we  have  studied 

7  8 

resonance  fluorescence  of  two-level  atoms  on  both  smooth  *  and  rough 

9  10 
surfaces  and  light  scattering  from  a  metal  grating. 

The  phenomenon  of  phase- conjugated  surfaces  (PCS)  has  been  recently 

explored,  where  both  the  propagation  direction  and  the  overall  phase  factor 

for  an  arbitrary  beam  of  light  may  be  precisely  reversed  via  nonlinear 
11  12 

optical  effects.  *  Such  effects,  in  general,  include  three-wave  mixing, 

e.g.,  second-harmonic  generation,  parametric  mixing  and  pockels  effects, 

(2) 

originating  from  the  x  term  of  the  susceptibility,  and  four-wave  mixing, 

e.g.,  third-harmonic  generation,  Raman  and  Brillouin  scattering,  dc  Kerr 

(3) 

effects  and  two-photon  absorption,  originating  from  the  x  terms  of  the 
susceptibility. 

In  the  present  paper,  we  shall  first  review  LISC  for  conventional 
surfaces  before  presenting  the  new  phenomena  on  PCS.  In  Section  2,  we 
review  resonance  fluorescence  of  two- level  atoms  near  smooth  and  rough 
surfaces.  A  radiative  dipole  near  a  PCS  will  be  studied  in  Section  3.  In 
Section  4,  laser- induced  processes  at  PCS  will  be  explored  in  detail,  where 
a  phenomenological  treatment  and  a  hydrodynamic  theory  for  laser- induced 
periodic  structure  and  PCS  reflectivity  are  analyzed.  Applications  of  LISC 
for  PCS  and  a  summary  are  presented  in  Section  5. 


II.  Resonance  Fluorescence  of  a  Two-Level  Atom  Near  a  Metal  Surface 


The  lifetime  of  an  excited  molecule  has  experimentally  been  found  to 

13 

vary  dramatically  as  a  function  of  distance  from  a  surface.  It  decreases 
or  increases  depending  on  the  distance  from  the  surface.  A  number  of 
researchers  have  examined  this  effect  from  the  viewpoint  of  reflected  field 
theory,  i.e.,  the  basic  calculation  is  concerned  with  the  interaction 
between  an  excited  molecule  and  its  own  reflected  radiation  field. 

Ref lected-f ield  theory  provides  generally  good  agreement  with 

.  ,  7-9,13-21 

experiments. 

When  an  adatom  is  driven  by  a  strong  resonant,  driving  coherent  field, 

it  creates  for  the  atom  or  molecule  an  environment  where  the  probability  of 

stimulated  emission  can  exceed  that  of  spontaneous  emission.  Under  this 

condition,  the  dynamic  ac  Stark  splitting  of  resonance  and  nutational 

oscillation  of  the  emitted  light  intensity  become  important  parts  of  the 

laser-driving  process,  such  that  interesting  "resonance  fluorescence"  and 

20 

other  nonlinear  optical  phenomena  can  occur. 

We  have  recently  derived  a  set  of  surface-dressed  optical  Bloch 

7-g  IQ-21 
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equations. 


by  which  we  can  examine  the  dynamics  and  relaxation  of 


an  adatom  on  the  surface.  In  these  equations,  we  are  able  to  evaluate  the 
resonance  fluorescence  spectrum  of  a  two- level  atom  near  a  metal  surface, 
taking  the  following  factors  into  account: 

(i)  The  reflected  electromagnetic  field  effects  due  to  the  presence 
of  the  interface 


(ii)  Collision  dephasing 

(iii)  Surface- induced  dephasing  due  to  the  reflected  photons 

(iv)  Resonance  excitation  of  surface  plasmons 

(v)  Random-phase  fluctuations  of  the  laser 


Based  on  these  studies,  we  have  also  been  extending  the  flat -surface 
resonance  fluorescence  calculations  to  the  rough  surface  case,  which  is 

9 

modeled  as  a  hemispheroidal  protusion  on  a  perfectly-conducting  surface. 

Ue  would  like  to  discuss  our  progress  on  the  surface  resonance  fluorescence 
spectrum  below. 


A.  Resonance  Fluorescence  at  a  Flat  Metal  Surface 

Recently  we  have  studied  the  resonance  fluorescence  spectrum  of  a  two- 

level  atom  near  a  metal  surface  by  means  of  a  set  of  surface-dressed  optical 

Bloch  equations  (SBE),  which  include  the  effect  of  surface-reflected 

photons,  i.e.,  photons  emitted  by  the  laser-driven  atom  and  reflected  by  the 
7  8 

surface.  '  Within  the  rotating-wave  approximation,  the  SBE  take  the  form 


s21(t> 

-Y2+iA 

in"(t)/2 

0 

§2i(t) 

0 

sr  «<t>  - 

m+(t) 

*Y1 

-ifl'(t) 

ft(t)  - 

1  ,  (1) 

S12(t) 

0 

-W+(t)/2 

-Y2-iA 

§12(12) 

0 

where  ft  =  |2><2| 

-  |  !><1 1 

in  the  population  inversion  of  the  atom. 

§..  is 

product  of  exp(ioj^t)  and  the  transition  operator  |i><j|,  £J±(t)  is  the  time- 
dependent  Rabi  frequency,  A  *  is  the  detuning  of  the  laser  (w^) 

with  respect  to  the  two-level  atom  (u^),  an<*  t*le  dephasing  rate  constant  y j 
is  the  sum  of  the  radiative  dephasing  y2  an<^  the  surface- induced  dephasing 
Ys  s  (2/!i)Im(f)|p12|2. 

To  calculate  y  ,  we  must  know  the  reflected  field  [Eq.  (1)}  or  the 
s 

complex  function  f(d),  which  can  be  determined  by  the  Sommerfeld-Hertz 
vector  procedure. ^  For  the  case  of  a  laser-excited  atom  emitting  near  a 
metallic  half-space,  the  imaginary  part  of  f(d)  has  the  form 


‘  '  n  (1  +  ej) 

x  {— [w2  -  ui2)2  +  P'l-  -  e2o)^] tnsinD  -  -^rcosD] 

E1  P  w2  1  D2 

+  2co6w2 [ r)cosD  +  — r  sinD]}  ,  (2) 

P  D2 

3 

where  C  *  1  and  n  5  ( 1  /D  )  +  (1/D)  for  the  case  of  the  induced  dipole 

3 

oriented  parallel  to  the  surface,  C  *  2  and  n  =  l/D  for  the  perpendicular 
case,  Ej  is  the  dielectric  constant  for  the  gas  medium,  k^  *  tu/e^/c  is  the 
wave  number,  u>  is  the  emission  frequency,  c  is  the  speed  of  light,  and  D  = 
2k^d  is  the  reduced  distance  which  is  dimensionless.  The  Drude  model'1  ^  has 
been  used  to  determine  the  dielectric  constant  of  the  metal  medium,  i.e.. 


e2U)  *  1  *  U(U  +  i«)  *  (3) 

where  6  is  the  inverse  of  the  relaxation  time  and  u  is  the  plasma 

P 

-1/2 

frequency.  We  note  that  Im(f)  has  a  resonance  at  u  *  ui  (1  +  e.)  , 

sp  p  l 

which  is  the  surface  plasmon  frequency  of  the  interface,  so  that  the 
reflected  field  given  by  Eq.  (1)  involves  resonance  energy  transfer  between 
the  adatom  and  the  excited  surface  plasmon. 

In  our  model,  the  driving  and  reflected  fields  have  been  treated 
semiclassically,  and  we  have  assumed  the  atom-surface  distance  to  be  large 
(>  30  nm).  The  relaxation  time  included  in  the  Drude  model  for  the  complex 
dielectric  constant  of  the  metal  medium  represents  the  dissipation  of 
electron  gas  which,  together  with  surface  plasmon  resonances,  influences  the 
surface  reflectivity  and  hence  the  behavior  of  the  reflected  field.  Here  we 
could  neglect  nonradiative  transfer  of  energy  from  the  excited  atom  to  the 


.*  .*  '  /  * 


metal.  Effects  of  the  laser  bandwidth  are  included  by  means  of  a  phase- 
diffusion  model  for  the  driving  field.  In  the  weak-field  or  large-detuning 
limit,  the  power  spectrum  of  scattered  light  has  two  peaks:  one 
corresponding  to  Rayleigh  scattering  at  the  laser  frequency  and  the  other 
to  fluorescence  at  the  atomic  transition  frequency  ^or  a  sufficiently 

strong  driving  laser  field,  the  spectrum  exhibits  three  peaks:  a  central 


one  at  (Rayleigh  component),  a  left  one  at  cu_  =  2oj^  -  -  6 


(three- photon  component)  and  a  right  one  at  w+  *  +  6  (fluorescence 

component),  where  6  is  the  ac  Stark  shift.  Results  for  the  peak  heights  H  , 

8 


Hq  and  H+  are  shown  in  Fig.  1  for  the  case  of  a  silver  surface.  The  key 
feature,  which  is  a  unique  behavior  due  to  the  surface,  is  that  for  certain 
atom-surface  distances  H_  is  larger  than  H+,  whereas  in  the  pure  gas-phase 
resonance  fluorescence  spectrum  H_  is  always  less  than  H+  (due  to  molecular 
collisions  for  positive  detuning).  Also,  the  population  inversion  of  the 
adatom  and  the  resonance  fluorescence  spectrum,  as  well  as  the  surface- 


induced  phase-decay  constant  of  the  adatom,  show  strong  oscillatory  behavior 

7,8 


as  a  function  of  the  adatom- surface  distance. 


B.  Resonance  Fluorescence  at  a  Rough  Surface 

For  a  rough  surface,  let  us  consider  a  surface  profusion  modeled  by  a 

9  22  23 

prolate  hemispheroid  on  top  of  a  plane.  '  ’  We  note  that  this  model  has 

24 

been  shown  to  be  identical  to  a  full  spheroid  in  a  vacuum,  so  that  our 
calculations  can  also  be  used  for  an  ellipsoidal  cluster.  The  two-level 
adatom,  which  is  located  at  a  distance  d  from  the  top  of  the  hemispheroid, 
is  driven  by  a  laser.  The  prolated  spheroidal  coordinates  (£,n»$)  are  used 
to  calculate  the  reflected  field.  We  define 


£  ■  (a  +  d)/f 


(4) 


«0  =  a/f 


(5) 


i—  *■«  *  «»  t 


f  »  (a2  -  b2)1/2 


where  a  and  b  are  the  semi-major  and  semi-minor  axes  of  the  hemispheroid. 

The  reflected  field  at  the  position  of  the  adatom  (transition  dipole)  in  the 

22 

near-field  approximation  can  be  written  as 


Er  '  ~t  + 


Mf^)' 


where  denotes  the  Legendre  function  of  the  second  kind  and  p  is  the 

induced  dipole  moment.  The  expansion  coefficient  depends  on  the  laser 

amplitude  and  parameters  and  Based  on  the  surface-dressed  optical 

7  8 

Bloch  equations  *  appropriate  for  the  excitation  and  dissipation  of  a  two- 
level  adatom  near  a  hemispheroid,  we  find  the  surface- induced  phase 
relaxation  constant  ys  to  be  given  by 

Y_  -  (2/h)Im(P)  ,  (8) 


where 


i  u-«>«0o;u.) 

l-r  LqjU,,)  -  «!«;(«„)  +  rl 


where  e  is  the  complex  dielectric  constant  and  r  is  a  somewhat  complicated 
function  of  Cq>  and  e. 

A  sharp  resonance  enhancement  of  the  adatom-hemispheroid  interaction, 
through  the  reflected  field  at  the  atomic  site,  occurs  when  the  specific 
shape  of  the  prolate  hemispheroid  corresponds  to  a  resonant  excitation  of 
plasmons  (see  Fig.  2).  The  time  oscillation  of  the  level  population 


■  *  -  - .  -  ■  -  V.  /.  t.  A  A  S  A  «*.  A  ••  • -  .  -  .  ■  -  .  •*.  •  . 
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decreases  as  the  shape  of  the  hemispheroid  approaches  the  plasmon  resonance. 
We  also  find  that  the  strong-field  three-peak  fluorescence  spectrum  is 
strongly  influenced  by  the  roughness  of  the  surface.  The  resonance 
excitation  of  the  plasmon  in  the  hemispheroid  remarkably  enhances  the 
adatom-surface  coupling,  such  that  the  dephasing  processes  broaden  the 
linevidths  of  the  spectrum.  In  the  small  detuning  case,  the  spectrum  has  a 
distinctive  three-peak  nature,  where  the  three-photon  side  peak  has  a 
measurable  height  which  is  almost  comparable  to  the  fluorescence  side  peak. 
The  plasmon  and  ref lected-f ield  broadening  influence  equally  the  three 
spectral  peaks.  In  the  large  detuning  case,  the  height  of  the  three-photon 
peak  is  decreased,  and  the  three-peak  spectrum  is  transformed  to  the  weak- 
field  two-peak  structure. 


III.  Radiative  Dipole  in  the  Presence  of  a  Phase-Conjugated  Surface 
A.  Phase-Conjugated  Mirror  Cavity 

Recently,  we  have  evaluated  the  lifetime  variation  of  a  dipole  centered 

26 

in  a  spherical  cavity  made  by  a  phase  conjugated  mirror  (PCM)  (see  Fig. 

3). 

27 

Following  Kuhn,  the  equation  of  motion  of  the  dipole  p  (assumed  to  be 
harmonically-bound  charge)  can  be  written  as 

2  2 

P  +  Wjj,__Er  -  bQ]u  ,  (10) 


where  u>  is  the  oscillation  frequency  in  this  absence  of  all  damping,  m  is 
the  effective  mass  of  the  dipole,  E^  is  the  reflected  field  at  the  dipole 
position  due  to  the  presence  of  the  PCM  cavity,  and  b^  is  the  damping 
constant  (inverse  lifetime)  in  the  absence  of  the  mirror.  The  dipole  moment 
p  and  the  (complex)  reflected  field  oscillate  at  the  same  frequency. 


-i[u>  +  Awjt  -bt/2 
M  *  Mq  e  e 


_  „0  -  i  [  u)  +  Aw]  t  -bt/2 

Br  =  Br  e 


where  and  Er  are  the  amplitudes  of  the  dipole  oscillator  and  reflected 

field,  respectively,  and  Aw  and  b  are  the  frequency  shift  and  the  inverse  of 

the  lifetime,  respectively,  in  the  presence  of  the  PCM  cavity.  The 

frequency  shift  is  found  to  be  quite  small  and  is  thus  unimportant  for  the 

2 

purpose  of  discussion  in  the  paper.  Recognizing  that  b  and  the  magnitude 

2  0  2 
of  (e  /pq1ii)Er  are  normally  very  small  compared  to  w  ,  we  then  have 


b  *  b_  +  [r-5—  ]  I  (ET) 
0  2u„mw  m  i* 


To  examine  the  variation  of  dipole  lifetime,  we  reduce  the  problem  to  the 

calculation  of  the  reflected  field  by  the  PCM  at  the  dipole  position. 

28 

By  using  the  Hertz  vector  method,  the  field  near  the  dipole  is  easily 
found  to  be 

E  -  2pnk*3  (n,  +  in,)[ - - — T - - — j]  COS0  e  lu>  t  (14) 

r  °  1  2  (k*r0)3  (k*rQ)2 

Eq  -  pQk*3  (o1  +  i02)I - - — 3 - - — sin6  e  ill)  1 

(k*rQ)  ^k*r0) 

where  Tq  is  the  linear  dimension  of  the  dipole  (normally  «  1),  and  n  * 

Hi  +  in2  the  complex  reflectivity  of  the  PCM  mirrors  where 


w*  ■  w  -  i2  =  ck*  , 


c  is  velocity  of  light,  and  the  definitions  of  E  ,  E  ,  6  and  r  are  readily 

r  u 

shown  in  Fig.  3.  We  note  that  b  <<  w,  and  taking  an  average  with  respect  to 


0  for  the  field  at  the  dipole  site,  we  obtain  the  imaginary  part  of  the 
reflected  field  as 


T  ,..0,  -  3  .  3r  n2  n2  , 

Im(V  -  pQk  [  -  ]  ,  (17) 

(krQ)  (krQ) 

where  k  is  the  propagation  constant,  k  =  ui/c.  The  normalized  decay  constant 
is 


2  2  “1 
where  |n|  =  +  n2  *s  t*le  of  the  reflectivity,  and  5  is  the 

normalized  lifetime  of  the  dipole. 

Figure  4  displays  typical  curves  for  the  variation  of  the  lifetime  for 
a  dipole  in  a  PCM  cavity,  as  a  function  of  with  a  fixed  reflectivity 

|  n  |  *  1%.  2Z  and  3%.  It  is  noted  that  due  to  the  focusing  nature  of  the 
reflected  radiation  wave  by  the  PCM,  the  coupling  between  the  dipole  and 
reflected  field  can  be  very  strong.  In  the  case  of  the  phase-shifted 
reflection,  i.e.,  n2  *  0>  there  are  three  cases  worthy  of  note: 

(a)  There  is  a  critical  value  of  n->/n,  *  defined  as  (n7/h,)  *  when  ru/h, 

l  l  l  l  c  11 

■*  (Hj/hj)  and  b  0,  i.e.,  in  this  case  the  lifetime  b  *  of  the  dipole  in 
the  PCM  cavity  can  be  infinite. 

(b)  In  the  region  n2/n^  >  0,  E  >  1,  the  dipole  in  the  PCM  cavity  can 
decay  faster  than  that  in  free  space. 

(c)  For  the  region  of  n2/n^  <  (n2/n^)c»  t^ie  dipole  decay  constant  E 


can  be  negative,  i.e.,  the  amplitude  of  the  dipole  oscillator  can  be 
amplified,  due  to  the  pumping  process  in  the  PCM. 


11 


IV.  Laser-Induced  Periodic  Structure  (LIPS)  on  Phase-Con lugated  Surfaces 

Laser-induced  periodic  structure  (LIPS)  on  conventional  surfaces  has 
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been  studied  over  the  past  several  years,  where  only  one  laser  was  needed 

for  LIPS,  which  has  been  explained  by  the  possible  mechanisms:  (i) 

formation  of  a  surface  polariton  with  periodic  electric  field  acting  on  the 
19  30  31 

surface  atoms;  ’  *  (ii)  laser-induced  longitudinal  acoustic  phonon  of 

32 

the  substrate  which  forms  a  standing  wave;  and  (iii)  small 

inhomogeneities  or  roughness  of  the  surface  layer  which  interact  with  the 

33 

incident  laser  and  forms  a  dipole  layer.  In  mechanism  (iii)  the 

interference  between  the  dipole  field  and  the  refracted  field  in  the 

substrate  in  turn  leads  to  inhomogeneous  energy  absorption  and  thus  the 

redistribution  of  the  surface  atoms. 

In  contrast  to  the  above  described  LIPS  on  a  conventional  surface, 

where  only  one  laser  beam  is  needed,  two  laser  beams  are  usually  required 

for  LIPS  on  PCS.  The  novel  method  of  LIPS  on  PCS  using  the  spatially 

nonuniform  reflectivity  variation,  both  in  amplitude  and  phase,  of  laser- 

influenced  interface  systems  include:  (i)  electron  concentration  changes 

(for  metal  and  semiconductor  surfaces);  (ii)  surface  deformation  due  to 

laser-induced  pressure  (or  density)  variation  of  the  surface  layer;  and 

(iii)  the  nonuniform  thermal  expansion  of  the  absorbing  layer.  Surface 

heating  and  deformation  by  a  single  pulsed  laser  have  been  reported  for  both 
34  35 

metal  and  semiconductor  surfaces.  Laser- induced  profiles  of  absorbing 

dye  in  solution  have  also  been  proposed,  where  two  interfering  laser  beams 

(at  the  same  frequencies)  were  used  for  LIPS  caused  by  the  concentration 

36 

modulation  of  the  dye. 

In  this  section,  we  »nall  present  a  unified  treatment  for  LIPS  on  PCS 
based  upon  a  hydrodynamic  theory.  The  mechanisms  of  LIPS  and  the  formation 


I 
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of  PCS  will  be  explored  mathematically.  Special  cases  with  analytic  results 
and  the  appropriate  conditions  for  the  validity  of  previously  published 
results  will  be  discussed. 

A.  Phase-Conjugated  Surface 

Phase  conjugation  has  been  conventionally  recognized  using  stimulated 
scattering  (Raman,  Brillouin  and  Rayleigh)  and  parametric  processes  (four- 
wave  mixing  and  two-wave  mixing). The  novel  method  of  phase  conjugation 
using  the  nonuniform  reflectivity  of  the  absorbing  layer  is  shown  in  Fig.  5, 
where  the  initially  uniform  reflecting  surface  may  form  a  periodic  structure 
caused  by  the  interference  of  two  incident  beams.  This  process  requires  the 
nonlinear  optical  effects  of  the  surface  layer,  e.g.,  nonlinear  polarization 
P^1  *  where  ^  is  the  total  wave  vector  given  by  the  sum  of  the 

reference  (strong)  wave  and  the  signal  (weak)  wave  by 

i(u>  t  +  1c. -r)  i(w_t  +  L-r) 

-  21  +  13  -  1  +  E3e  +  c.c.)  .  (19) 

The  laser- inf luenced  reflectivity  (R),  which  in  general  is  time-dependent 
and  is  a  complex  variable  with  amplitude  and  phase  modulation,  may  be 
expressed  (up  to  the  first-order  expansion  of  the  laser  intensity)  as 

R  -  Rq  +  PHj  ,  (20) 

where  Rq  is  the  reflectivity  from  a  conventional  surface  and  is  proportional 
to  the  total  intensity,  |E^|  +  | |  .  p  is  the  expansion  coefficient  of 

the  interference  pattern  (or  reflectivity  modulation)  formed  by  the  incident 


waves,  and  R^  is  given  by 


where  we  have  assumed  the  degenerate  case,  ■  w  ,  and  have  used  the 
rotating-wave  approximation  (RWA).  Ak  is  the  phase  mismatching  vector 
defined  by  Ale  =  •  For  the  one-dimensional  and  degenerate  (k^  *  * 

k)  case  |A^*r|  =  2kzsin(8/2),  such  that  the  interference  pattern  is  governed 
by  the  angle  (6)  between  the  incident  waves.  We  shall  show  in  the  next 
section  that  the  lifetime  of  the  interference  pattern  is  proportional  to  the 
inverse  of  |Alc-r|2,  i.e.,  a  small  angle  (in  the  order  of  few  mrad)  is 
required  for  a  long-lived  interference  pattern  to  be  observed. 

The  overall  reflected  waves  defined  by  the  laser-modulated  reflectivity 
may  now  be  calculated  by  using  RWA  to  eliminate  the  irrelevant  terms, 
whereby  we  obtain  from  Eqs.  (19)  and  (20) 

Kef  '  +  *3> 
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Thus  the  total  reflection  field,  f  consists  of  three  parts:  (i)  the 
usual  reflected  fields  from  the  the  unperturbed  surface,  the  Rq  terms;  (ii) 
the  reflected  fields  from  the  laser-perturbed  surface  in  the  usual  way, 
i.e.,  from  the  incident  r  to  the  reflected  -r  direction;  and  (iii)  the 
reflected  fields  in  the  phase -conjugated  forms,  the  E*  (r)  terms.  The 
reflectivity  of  the  phase-conjugated  surface  (PCS)  defined  by  the  intensity 
ratio  between  the  signal  wave  and  its  conjugation  becomes 

RpCS  -  P^EJ4  -  (SI^2  .  (23) 

This  is  proportional  to  the  square  of  the  nonlinear  coefficient  of  the 
surface  layer  and  the  intensity  square  of  the  reference  wave  (E^),  which  is 


X  s 


usually  much  stronger  than  the  signal  wave  (E^).  We  note  that  the  above 
obtained  results  involved  several  assumptions:  (i)  the  pump  depletion  of 
the  reference  wave  is  ignored;  (ii)  steady-state  reflectivity  (or 
interference  pattern)  is  assumed;  (iii)  instant  energy  transfer  (or  local 
heating)  of  the  surface  layer  is  assumed;  and  (iv)  heat  diffusion  of  the 


S&S 


m 
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surface  layer  and  attenuation  of  the  fields  are  ignored.  For  a  rigorous 
treatment,  these  assumptions  will  be  removed  in  the  following  sections. 

B.  Hydrodynamic  Theory  for  LIPS 

The  generation  of  density  and  temperatures  fluctuations  following  the 


absorption  of  light  is  known  as  a  stimulated  thermal  process,  where  the 

light-induced  modulation  of  the  refraction  index  of  the  absorbing  layer  can 
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act  as  a  phase  grating  in  the  diffration  of  incident  light.  The  full 
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dynamics  of  a  thermal  grating  may  be  described  by  a  hydrodynamic  theory 

which  includes  the  thermal  processes  of  the  absorbing  media  and  the  optical 

nonlinear  effects  (or  polarization)  caused  by  laser  fields.  The  complete 

set  of  hydrodynamic  equations  which  couple  the  thermodynamic  variable 

(density,  internal  energy  and  temperature)  of  the  medium  and  the  electric 
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fields  of  the  incident  beams  are  given  by 

-  —  I2  +  DV2!)  -  U/t  ,  (24) 

°v  If  -  k,2t  -  Fff  ■  u/t  •  (25) 
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Here  U  Is  the  Internal  energy  of  the  medium  with  absorption  coefficient  a, 

diffusion  constant  D  and  thermalization  time  x;  T  and  p  describe  the 

temperature  and  density  fluctuation  (normalized  by  their  equilibrium 

values);  K  is  the  thermal  conductivity  and  is  the  specific  heat  at 

constant  volume;  F  is  an  appropriate  coupling  constant  and  q  is  the 

effective  viscosity;  B^,  is  the  bulk  modulus  (or  thermal  expansion 

coefficient);  and  y  is  the  electrostrictive  constant.  Equation  (26)  is 
e 

just  the  usual  Maxwell  equation  with  a  nonlinear  polarization  relating  to 
the  dielectric  constant  (c),  density,  temperature  and  field  amplitudes  by 


-  [(f^)T  P  +  (|f)p  T](S/4it)  =  (Yep  +  yTT)(f/4ir)  .  (28) 


Employing  the  slowly-varying  envelop  approximation  (SVEA)  and  RWA  which 
eliminates  the  irrelevant  (or  nonresonant)  terms,  we  obtain  a  set  of 
simplified  equations  for  the  amplitudes  of  the  density,  temperature  and 
electric  fields  as 
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(29) 

(30) 

(31) 

(32) 

(33) 


where 


are  the  k-vector  and  decay  time  of  thermal  grating  and  the  specific  heat 

ration,  respectively;  and  are  the  frequency  and  linewidth  of  the 

2 

acoustic  phonon  wave  given  by  =  wave  velocity  *  q,  Tq  ■  nq  •  In  deriving 
the  above  equations,  we  have  further  assumed  the  degenerate  case,  * 

u,  kj  ■  kj  *  k  ■  nw/c  and  ignored  the  small  absorption  in  the  thin  layer, 
exp(a)  =  1.  The  amplitudes  of  the  incident  (t^),  reflected  (^)*  probing 
(&j)  and  the  phase -conjugated  (1^)  fields,  internal  energy,  density  and  the 
temperature  are  defined  by 

i(w. t  -  it  r) 

J  +  c.c.)  ,  (35) 

where  i))^  =  E . ,  p^,  T^,  with  j  *  1,2, 3, A.  Similar  to  that  of  and  E^, 
equations  for  E^  and  E2  are  not  shown  here.  They  are  needed  only  when  the 
pump  depletion  i.e.,  the  variation  of  E^  with  z,  is  considered. 

We  note  that  within  the  SVEA,  the  thermal  grating  is  formed  in  such  a 
way  that  only  the  relevant  component  of  the  reflected  wave  meeting  the 
energy  and  momentum  conservation  is  strongly  coupled,  where  the  phase 
matching  condition  is  met  between  the  the  probe  (E^)  and  the  conjugated  (E^) 
waves  through  this  intrinsic  feature  of  the  thermal-grating-induced  acoustic 
wave.  This  feature  may  be  clearly  realized  by  the  steady-state  solutions 
for  the  coupled  equations  ( 29 ) - ( 33)  which  reduce  to 


5 r  '  -’ll'll E3  -  ®2E1E2EJ 


2 r  ■  %  +  Wi'3 


where  are  the  appropriate  coupling  constants,  proportional  to  the  decay 
time  of  the  thermal  grating  defined  in  Eq.  (3A).  These  steady-state 
equations  are  similar  to  those  of  the  degenerate  four-wave  mixing  FWM), 


where  E^  and  E^  are  the  pump  beams  and  E^  and  E^  are  the  signal  (probe)  and 
the  conjugate  beams.  The  first  terms  on  the  right-hand  sides  of  the  above 
equations  describe  the  nonlinear  refractive  index  changes  caused  by  the 
optical  Kerr  effects,  while  the  second  terms  account  for  the  phase- 
conjugated  FV)M.  A  complete  description  of  the  dynamics  of  LIPS  requires  the 
numerical  solutions  of  the  above  described  coupled  equations.  For  tractable 
results,  we  shall  discuss  some  of  the  special  cases  under  appropriate 
assumptions. 

C.  Special  Cases 
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For  the  linearized  case  (or  the  liquid  limit)  ,  the  absorbed  energy  of 
the  surface  layer  is  instantly  and  locally  thermalized,  i.e.,  the  diffusion 
term  is  neglected  and  the  internal  energy  reaches  its  steady-state  obtained 
by  taking  the  limit  of  D  +  0  and  t  >  0  in  Eq.  (29), 


U  =*  It2£(E.E*  +  E*E. )  . 

s  4  n  1  3  2  4 


In  this  case  Eq.  (30)  becomes  a  usual  heat  diffusion  equation 


(at  +  rT  '  D“2)T1  "  Fl^l  "  S  ’ 
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where  we  have  generalized  the  heat  diffusion  equation  into  the  three- 

dimensional  system  where  the  heat  diffusion  in  the  radial  direction  causes 

2 

the  decay  of  the  thermal  grating  given  by  *  Dq  ,  and  heat  absorption  of 
the  surface  layer  in  the  normal  (z)  direction  causes  the  LIPS  via  the  non- 
uniform  thermal  expansion  of  the  surface  layer.  F^  is  a  coupling  constant 
and  denotes  for  the  time  derivative  of  the  density.  We  shall  now 
calculate  the  amplitude  (height)  of  LIPS  based  on  a  simplified  version  of 


Eq.  (39)  called  the  "surface  heating"  system.  This  reduces  the  "volume 

n  / 

heating"  system  of  Eqs.  (39)  and  (40)  intoJ 


(3t  +  rT)Tl  "  D7 ~2  T1  +  Fl^l  ' 
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(41) 


with  the  initial  condition  T^(t  ■  0)  =  0  and  boundary  condition 


3T 


J 1.  -  -(1  -  R)En  , 

3z  >z  =  0 


(42) 


where  R  is  the  usual  reflectivity  of  the  incident  E^  ^  and  the  reflected 

E_  .  wave,  and  En  is  defined  by  E_  s  E.E*  +  EiE.  evaluated  on  the  surface 
3,4  0  •'01324 

plan,  z  =  0.  Under  the  condition  of  isotropic  thermal  expansion,  the 
amplitude  of  LIPS  is  given  by 

Q(t)  =  y  j”dz  T^z.t)  .  (43) 

Thus  the  variation  of  Q(t)  is  governed  by 

(§£  +  rT)Q(t)  *  Bt(1-R)Eq/3  ,  (44) 

where  we  have  neglected  the  density  coupling  term,  and  in  deriving  Eq.  (44), 
we  have  used  Eqs.  (41)  and  (42)  and  also  assumed  (3T^/3z)  at  z  »  »  is  zero. 
The  amplitude  (height)  of  LIPS  may  then  be  calculated  from  Eq.  (44)  whose 
exact  solution,  for  Q(T  ■  0)  *  0,  is  given  by 

( l_R)Bm  rt  T  (t’-t) 

Q(t)  -  — 3 - -  J^dt*  e  T  EQ(t')  .  (45) 

We  note  that  the  time-dependent  profile  of  the  LIPS  depends  upon  the  thermal 
grating  decay  rate,  and  the  evolution  of  the  fields  Eg(t)  =  E^(t)E*(t)  + 

E^(t)E^(t)  evaluated  at  the  interface  plan  z  =  0.  For  the  cw  or  square 
pulse  case,  it  is  easy  to  find 


rT  -  Dq^  *  D(^) sin  (|)  ,  (46b) 

where  the  steady-state  amplitude,  for  T^t  >>  1,  is  proportional  to  E^/r^. 

We  note  that  formation  of  LIPS  with  significant  amplitude  requires  two 

conditions:  high  enough  applied  field  (Eq)  and  long  enough  thermal  grating 
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lifetime,  i.e.,  small  damping  factor  T^,.  Because  is  proportional  to  Dq  , 
we  require  a  small  diffusion  coefficient  (D)  and  a  small  angle,  i.e.,  q  * 
2ksin(^)  is  small.  The  latter  condition  may  be  achieved  by  lasers  with 
short  wavelength  and/or  by  small  propagation  angle  between  the  pump  (^)  and 
the  probe  (2^)  field. 

To  include  the  density  fluctuation  effects  in  LIPS,  the  following 
coupled  equations  may  be  obtained 


P  +  rpP  +  OqP  +  F2Q  -  BjEq  , 


Q  +  rTQ  +  FjP  *  b2e0  , 


(47) 

(48) 


wherewhere  P  is  the  density  fluctuation  (or  laser-pressure)  induced 
amplitude  of  LIPS,  and  are  the  associated  damping  factor  and 
frequency,  and  F^  2  and  B^  2  are  the  appropriate  coupling  constants.  The 
steady-state  solution  of  Eqs.  (47)  and  (48)  is  straightforward  and  will  not 
be  shown  here.  We,  however,  note  that  the  overall  amplitude  of  LIPS  in 
general  is  a  complex  variable  and  time -dependent. 

D.  Reflectivity  of  PCS 

2  2 

The  reflectivity  of  PCS  defined  by  Rp^g  =  |E^|  /|E^|  evaluated  on  the 
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z  *  0  plane,  in  general,  is  time -dependent.  Again,  numerical  methods  are 
required  for  the  transient  values.  For  the  case  of  a  strong  pump,  with  E.  _ 


»  E  f  the  steady-state  equations  (35)  and  (36)  for  the  case  of  p.  -  << 

1 1  3 

e2  A  and  P2  =  e4  -  g  yield 
E.(0>  ,  , 

Rpcs  =  Ie^ToT*  =  ta" "(|PEi(0)E2(°)l|)  ,  (49) 


where  p  is  proportional  to  the  decay  time  of  the  thermal  grating  defined  in 
Eq.  (34)  and  L  is  the  effective  thickness  of  the  absorbing  layer  for  the 
phase-conjugated  back  scattering  to  occur.  We  note  that  in  the  low- 
reflectivity  limit  and  E^(0)  =  E^(0)  =  1^,  Eq.  (49)  reduces  to  Eq.  (23)  when 
L  *  1.  Furthermore,  the  PCS  reflectivity  can  actually  be  much  higher  than 
one,  i.e.,  the  phase-conjugated  wave  (E^)  is  coupled  to  the  strong  pump 
fields  (E^  and  E2)  and  its  intensity  is  "amplified"  to  a  level  much  higher 
than  that  of  the  probe  field  (E^).  The  infinite  peaks  are  predicted  from 
Eq.  (49)  when  the  conditions  of  tan  x»(n+J)ir  are  met.*1  We  should, 
however,  note  that  this  behavior  occurs  only  at  the  weak  signal  limit  and 
does  not  occur  when  the  pump  depletion  is  included  (strong  signal  case). 

5.  Applications  and  Summary 

In  this  paper  we  have  reviewed  the  resonance  fluorescence  of  two- level 
atoms  on  both  smooth  and  rough  surfaces.  For  a  system  of  smooth  surface,  we 
have  shown  the  oscillatory  behavior  for  both  the  resonance  fluorescence  and 
the  adatom  lifetime  as  a  function  of  the  adatom-surface  distance  (referred 
to  Fig.  1).  For  an  adatom  on  a  rough  surface,  we  found  that  the  sharp 
resonance  enhancement  caused  by  the  plasmons  excitation  occurs  for  a 
specific  shape  of  the  prolate  hemispheroid  (see  Fig.  2).  Results  obtained 
from  this  system  can  provide  us  useful  characterization  in  the  area  of 
surface  chemistry,  such  as  the  orientation  and  structure  of  the  adspecies, 
the  strength  of  the  adatom-surface  coupling  and  the  role  of  surface 
roughness  on  the  lifetime  and  fluorescence  of  the  adspecies. 
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For  phase -conjugated  systems,  we  have  proposed  a  new  phenomenon  in 

Section  3,  of  dipole  radiation  on  a  PCS  where  the  dipole  lifetime  can 

actually  be  infinite.  In  a  pumping  process,  the  reflected  field  may  be 

"amplified"  by  a  set  of  counter-propagating  strong  fields.  In  Section  4,  we 

have  discussed  LIPS  caused  by  the  interference  between  two  incident  waves. 

The  phase -conjugated  field,  originated  from  the  optical  nonlinear  effects  of 

the  absorbing  layer,  may  be  achieved  with  a  reflectivity  much  larger  than 

one.  Applications  of  LIPS  include,  for  example,  (i)  formation  of 

microstructures  for  microelectronic  fabrication,  where  the  spacing  between 

the  periodic  gratings  is  only  limited  by  the  laser  wavelength  and  the 

incident  angles  of  the  beams  [Eq.  (46b)];  (ii)  for  the  application  of  real 

time  holography,  greater  resolving  power  is  necessary,  i.e.,  greater  PCS 

reflectivity  operated  at  a  larger  angle  (6).  This  may  be  achieved  by  higher 

incident  beam  intensities  and/or  shorter  wavelengths.  The  materials 

suitable  for  the  above  described  applications  will  be  an  absorbing  dye  with 

very  high  absorptibility  and  very  low  thermal  conductivity  but  high  bulk 

expansion  coefficient,  referred  to  Eq.  (45).  An  experimental  study  of  LIPS 

using  a  He-Ne  laser  and  an  absorbing  layer  of  epoxy  resin  mixing  with  a 

40 

brilliant,  green  dye  has  been  reported. 
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Figure  Captions 

1.  Atom-surface  distance  dependence  of  the  heights  of  the  three  peaks  in 

the  resonance  fluorescence  spectrum  for  the  case  of  a  silver  surface, 
where  H  ,  and  H+  correspond  to  the  peak  heights  (in  arbitrary  units) 
of  the  three-photon,  incoherent  Rayleigh  and  fluorescence  components, 
respectively.  D  =  2kd  is  the  reduced  distance,  where  k  is  the 
wavenumber  and  d  is  the  atom-surface  distance.  The  Drude  model  is  used 
for  the  dielectric  function  of  the  metal,  while  the  dielectric  function 
of  the  gas  medium  is  set  equal  to  one.  The  detuning  is  1  A  and 

the  Rabi  frequency  is  10  A,  in  the  unit  of  Einstein's  A-coef f icient  for 
spontaneous  decay.  The  solid  curves  are  for  the  induced  transition 
dipole  of  the  atom  oriented  perpendicular  to  the  surface,  and  the  dashed 
curves  are  for  the  parallel  case. 

2.  Rough  surface- induced  phase-decay  constant  y  as  a  function  of  the  semi- 

s 

minor  axis  b  of  the  surface  protrusion,  where  the  semi-major  axis  is 

fixed  at  100  A.  y  is  in  the  unit  of  Einstein's  A-coeff icient. 
s 

3.  Geometry  of  an  electric  dipole  in  a  phase-conjugated  (time-reversal) 
cavity.  The  dipole  with  moment  y  =  ye^  is  located  at  the  center  of  the 
cavity,  where  e^  is  the  unit  vector  along  the  x^-axis.  The  emitted 
radiation  field  at  the  point  (r,9)  is  denoted  by  (E  ,E.)  inside  the 
cavity. 

4.  Normalized  decay  rate  constant  5  of  the  dipole  versus  with  fixed 

reflectivity  |q|:  (a)  |n|  =  1Z;  (b)  |n|  *  2Z  (c)  |n|  =  3Z.  The  results 

show  that  the  lifetime  of  the  dipole  is  sensitive  to  |q|  and  r^/n^  of 
the  PCM  cavity.  There  is  a  critical  value  of  that  leads  to  5  ■  0, 

i.e.,  the  lifetime  approaches  ii finity.  In  these  curves,  the  values  of 


A  =  0.3  ym  and  Tq  =  100  A  are  used. 
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